Electronic thermal conductivity at high temperatures: Violation of the 
Wiedemann-Franz law in narrow band metals 
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We study the electronic part of the thermal conductivity k of metals. We present two methods 
for calculating «, a quantum Monte-Carlo (QMC) method and a method where the phonons but 
not the electrons are treated semiclassically (SC). We compare the two methods for a model of 
alkali-doped C60, A3C60, and show that they agree well. We then mainly use the SC method, which 
is simpler and easier to interpret. We perform SC calculations for Nb for large temperatures T 
and find that k increases with T as ft(T) = a + bT, where a and b are constants, consistent with a 
saturation of the mean free path, I, and in good agreement with experiment. In contrast, we find 
that for A3C60, K (T) decreases with T for very large T. We discuss the reason for this qualitatively 
in the limit of large T. We give a quantum-mechanical explanation of the saturation of I for Nb and 
derive the Wiedemann-Franz law in the limit of T <C W , where W is the band width. In contrast, 
due to the small W of A3C60, the assumption T < W can be violated. We show that this leads to 
k(T) ~ T- 3/2 for very large T and a strong violation of the Wiedemann-Franz law. 
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I. INTRODUCTION 

In the Sommerfeld theory, the electronic contribution 
k to the thermal conductivity of a metal can be written 



1 

K = -v F lc v = — 

6 6 mvF 



— 



(1) 



where vf is the Fermi velocity, I the mean free path, 
c v = tt 2 kgnT ' I (mvp) is the specific heat per electron, n 
is the electron density, ks is the Boltzmann constant, m 
is the electron mass and T the temperature. For large T 
the mean free path decreases as I ~ and one might 
then expect that k approximately approaches a constant 
for large T. This is indeed found for many good metals. 
However, in the context of electrical conductivity, it has 
been found that for large T the so-called parallel resistor 
formula describes the experimental resistivity quite well 
for many metals with a large resistivity— This formula 
corresponds to the assumption 



I = d+ —. 
T 



(2) 



where d is a distance of the order of the separation of 
two atoms and c is a constant. According to this formula 
the mean free path initially decreases rapidly with T but 
then saturates at I ~ d. This formula is often justified 
by arguing that in a semiclassical picture, at worst, an 
electron is scattered at every atom, leading to I ~ d. 
Inserting Eq. (0) in Eq. (JIJ leads to 



a + bT, 



(3) 



where a and b are constants. This formula gives a fairly 
good qualitative description of experimental results for 
many transition metal and transition metal compounds 
with a small k. 

The theoretical justification for Eq. (0), however, is 
very unsatisfactory. It is based on a semiclassical pic- 
ture, which is only valid for I 3> d, and which cannot be 



used to discuss what happens for I ~ d. Nevertheless, in 
the 1970's and early 1980's it seemed that I > d was a 
universal behavior, 3 satisfying the Ioffe-Regel condition. 4 
Later experimental work, however, has found many ex- 
amples of metals where the resistivity is much larger than 
predicted by the Ioffe-Regel condition and where the ap- 
parent mean free path is much shorter than the separa- 
tion of two atoms£ This illustrates that the semiclassical 
explanation for Eq. is incorrect and there is a need 
for a quantum- mechanical explanation of Eq. Q. Such 
a theory has been presented in the context of the electri- 
cal conductivity^ The purpose of this paper is to provide 
a quantum-mechanical justification for the thermal con- 
ductivity in Eq. @- 

To calculate the thermal conductivity we use a Quan- 
tum Monte-Carlo (QMC) method, which can solve the 
models used here accurately. To interpret the results, 
we also introduce a simpler semiclassical (SC) method, 
where the phonons but not the electrons are treated semi- 
classically. 

For many metals, the Wiedemann-Franz law 



K _ 7T 2 ,fcg 2 



(4) 



is approximately satisfied, where a is the electrical con- 
ductivity and e the electron charge. We derive this law 
by making assumptions which are expected to be reason- 
able when T is so large that the transport is completely 
incoherent but T is still much smaller than W, where W 
is the band width. For T > W, however, we find that 
the Wiedemann-Franz law is strongly violated. This may 
apply to alkali-doped fullerides, A3C60, at or somewhat 
above the highest temperatures that can be achieved ex- 
perimentally. 

In Sec. [H] we present the models and in Sec. IIIII the 
methods used. In Sec. IIVI we compare the QMC and 
SC methods. The results for Nb and A3C60 in the SC 
method are given in Sec. These results are discussed 
qualitatively in Sec. I VII 
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II. MODELS 

We first consider a model of Nb, referred to as the 
transition metal (TM) model, which is appropriate for de- 
scribing transition metals or transition metal compounds. 
Each Nb atom has a five-fold degenerate (Nd = 5) level. 
The hopping matrix elements are described by where 
v = (m, i) is a combined label for a orbital index m and 
a site index i. Thus the electronic Hamiltonian is 



fl<7 



,t 



(5) 



where cj^ creates an electron in the state \i with spin 
a. We assume that the Coulomb interaction can be ne- 
glected, since the band width of Nb is rather large. The 
precise form of the hopping matrix elements has been 
described elsewhere^ 

We assume that the electron scattering is due to the 
electron-phonon coupling. For the Nb model, we assume 
that the phonons couple to hopping integrals (HI). The 
phonons are approximated as Einstein phonons, with one 
phonon for each coordinate direction. The frequency 
uiph = 0.014 eV was set equal to the average phonon fre- 
quency of Nb metaliS Due to the vibrations of the atoms 
the hopping matrix elements are modulated, leading to 
an electron-phonon coupling^ 

We next introduce a model of alkali-doped Ceo, A3C60, 
referred to as the Cgo model. We use a model includ- 
ing the partly occupied three-fold degenerate t\ u orbital 
on each Ceo and the hopping matrix elements connect- 
ing these orbitals. This results in a Hamiltonian of the 
same general form as in Eq. (JSJ but with a different lat- 
tice structure and different orbital degeneracy (N ( i = 3). 
The hopping integrals are obtained from a tight-binding 
descriptionpSai For simplicity, the Coulomb interaction is 
neglected, although this may be a questionable approxi- 
mation. The main electron-phonon coupling is due to the 
intramolecular five-fold degenerate phonons of H g sym- 
metry, which have an on-site Jahn- Teller coupling to the 
t lu levels of a Cgo moleculeiS We refer to this as a level 
energy (LE) coupling. For the TM and Ceo models we 
define the dimensionless coupling A so that the real part 
of the electron-phonon part of the electron self-energy is 
given by 



ReS cp (cj) = — Aoj, 



(6) 



in the weak-coupling limit and for small w. 

To introduce current operators we follow Mahan^fi 
The particle current is obtained from particle conserva- 
tion 

j = Ti [H ' P] = " \ ^ (Rm " R -)V^cw, (7) 

where H is the full Hamiltonian and P is the polarization 
operator 



(8) 



Here R„ = R; is the position of the ith atom (v = (m, i)). 
In a similar way the energy current is obtained from en- 
ergy conservation 



j B = ~^Ri[ffA] 



(9) 



where hi is a Hamiltonian of the ith atomic site, defined 
in such a way that the hopping terms between two sites 



are split equallybetween these two sites and H — ^ hi. 

(10) 



iyb< 

Following Ref. [l(J, we also introduce a heat current 



JQ = JE - AJ, 
where /1 is the chemical potential. 

III. METHODS 



The conductivity is calculated using a Kubo formalism. 
We assume an isotropic system and define a heat current 
- heat current correlation function 10 



r (22) 



(iuj n ) 



1 



3Nfl 



dr e -» T (r T j Q (T)-j Q (0)), (11) 



where N is the number of atoms, £7 is the volume per 
atom, (3 — l/(fcsT), T T is a time-ordering operator and 
uj n is a Matsubara frequency. In a similar way we define a 
particle current - heat current correlation function ■k^ 12 ^ 
and a particle current - particle current correlation func- 
tion 7r( n ). These correlation functions are analytically 
continued, giving retarded response functions tt™^(lu). 
Then the electrical conductivity is given byifi 



The heat conductivity i; 



„2 Im7r r ( et' 



huj 



,10 



Thuj 



r (22) 



H- 



[Iimr^V)] 2 

Im^ret 3 H 



(12) 



}. (13) 



We use two different methods for performing calcu- 
lations, a determinantal quantum Monte-Carlo (QMC) 
method—^ and a semiclassical (SC) method. The QMC 
method can be used to obtain properties very accurately. 
To interprete these results we use a SC method where the 
phonons but not the electrons are treated semiclassically. 

Due to the lack of a repulsive Coulomb interaction in 
the models treated here, the QMC method has no so- 
called sign problem. The statistical errors in the QMC 
calculation of response functions for imaginary times can 
then be made arbitrarily small by improving the sam- 
pling. To analytically continue these response functions 
to the real frequency axis we use a maximum entropy 
method^ 
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The Hamiltonians treated here can be written in the 
form 

fivna n 

= H cl (q)+H B , (14) 

where g^ un is a coupling constant, q n is a phonon coordi- 
nate, p n a phonon momentum and w p h is the frequency of 
the phonons. In the QMC calculation the starting point 
is the partition function 

Z = Tr|dq(q| e -^|q) (15) 

where Tr is a trace over all electronic states and 
the integral over q = (qi, q%, ...) results in a trace 
over q. Since H cl (q) and H B do not commute, we 
use a Trotter decomposition and define At — (3Ti/L, 
where L is some integer. We write exp(—f3H) as 
n^ 1 exp(— AriJ e ;(q))exp(— AtH b ). By inserting com- 
plete sets of phonon states between each factor, one 
obtainsi±*±£ 



where j% = j, j% — Jq and /; is the Fermi function for 
the energy We have here for simplicity considered 
the xx component. The result is averaged over different 
configurations of random displacements of the atoms. 

The electron-phonon interaction leads to a renormal- 
ization of the phonon frequency. This renormalization 
acts back on the electronic properties. This is included in 
the QMC approach above, but not in the SC method. To 
be able to compare more directly with this SC method, 
we also perform QMC calculations where the electron- 
phonon interaction is neglected when calculating the 
weight functions entering the sampling in Eq. (|16f) . In all 
other parts of the QMC calculation the electron-phonon 
interaction is fully taken into account. The result is that 
the phonon coordinates are sampled as if the phonons 
were free phonons. In this QMC calculation, the phonons 
are then not renormalized, and the corresponding influ- 
ence on the electronic properties is absent. 



IV. COMPARISON OF QMC AND 
SEMICLASSICAL METHODS 



Z = Tr y(nf =1 d qi )e- A ^^)( qi 



e- ArHB \q L ) 



.e- ATH ^(q 2 \e 



AtH b 



qi 



(16) 



The phonon coordinates q^ = {q n ,i} = {ln(ri)} are sam- 
pled by using a Monte-Carlo approach with the inte- 
grand as a weight factor. Here = iAt is an imaginary 
time. Correlation functions (X(Ti)X{rj)) are obtained 
in a similar way by inserting the operators at positions 
corresponding to r t and TjU'H For the case of the HI 
coupling there are technical complications in terms of 
updating the relevant quantities after each Monte-Carlo 
step, which are discussed elsewhere^ 

In the SC treatment of the phonons^ we introduce 
a super cell with periodic boundary conditions. Each 
phonon coordinate is given a random displacement ac- 
cording to a Gaussian distribution centered at zero and 
with the width 



Mui ph 



[n B (T) + 0.5] 



where 



n B (T) 



1 



Jluj ph /(k B T) _ ^ ' 



(17) 



(18) 



is the occupation of the phonon mode. For fixed phonon 
coordinates, the resulting Hamiltonian is a one-particle 
Hamiltonian, which can easily be diagonalizcd. The re- 
sponse functions 7r(" m ) can then be expressed in terms of 
the eigenstates \l) and eigenvalues Ei of this Hamiltonian 



ImTrir'M = -^Y(l\(jmU')(l\(3nU) 
ll' 



x(fi- fi>)5{huj + si -sj), 



(19) 



To test the accuracy of the SC method, we compare 
the thermal conductivity of the Cgo model according to 
the QMC and SC method. As discussed in Sec. IIIII 
in the QMC calculation the phonon frequency is renor- 
malized by the electron-phonon interaction, while this 
is not the case in the SC calculation. To compare the 
two methods we use two different approaches. In one 
approach we suppress the phonon renormalization in the 
QMC calculation, as discussed in Sec. IIIII In a second 
approach we calculate an effective phonon frequency, w G ff , 
in the QMC program and then put in this frequency by 
hand in the SC calculation, keeping the coupling con- 
stants fixed. To obtain the effective phonon frequency 
we calculate the phonon Green's function D(t) for imag- 
inary times r. We compare this with the phonon Green's 
function, Dq(t, u> e s), for noninteracting phonons with the 
frequency o; e ff- We minimize 



(20) 



with respect to w c ff, where Tj correspond to the discrete 
values of r used in the QMC calculation. The electron- 
phonon coupling constants are kept unchanged. 

In the SC method k(oj) is calculated directly for real 
frequencies, while in the QMC method we calculate the 
response functions 7r" m (r) for imaginary times and then 
use a maximum entropy method (MEM) to perform 
the analytical continuation. To make the QMC and 
SC methods more comparable, we have therefore in the 
SC method also made a transformation of the data to 
imaginary times for each random configuration, which is 
an accurate and well-behaved transformation. We have 
then transformed the data back to real frequencies, us- 
ing MEM. These results, referred to as SCMEM results, 
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FIG. 1: a) Thermal conductivity k(oj) of A3C60 as a func- 
tion of lu according to QMC (full line), SC (dashed line) and 
SCMEM (dashed-dotted line), b) The same as in a) but with 
the renormalization of the phonon frequency suppressed. The 
parameters are io p h = 0.1 eV, A = 0.5 and T = 0.1 eV. In a) 
uiph is renormalized to w e g = 0.053 eV. The model of A3C60 
described in Sec. HTIwas used. The calculations were done for 
a cluster with 3x4x4 sites and averaged over 20000 different 
random configurations. A Gaussian broadening with FWHM 
0.03 eV was used. The figure illustrates that the QMC and 
SCMEM calculations agree very well and that the downturn 
at lo — in the SC results is lost in the SCMEM results due 
to the MEM procedure. 



then contain errors introduced by the analytical continu- 
ation, and are in this sense more comparable to the QMC 
results. 

Fig. H shows results for k(oj) of the Cgo model for 
T = 0.1 eV, Lu ph = 0.1 eV and A = 0.5. In the 
QMC calculation the phonon frequency is renormalized 
to uj c s = 0.053 eV. In Fig. the phonon renormal- 
ization is taken into account and in Fig. it is sup- 
pressed. It is interesting that the QMC and SCMEM 
data agree quite well, both with and without renormal- 
ization of the phonon frequency. This suggests that the 
SC method is rather accurate. We also notice that the 
renormalization of the phonon frequency leads to an ef- 
fectively stronger electron-phonon coupling and therefore 



a smaller k(uj = 0). This follows since the coupling goes 
as the inverse phonon frequency if the coupling constants 
are kept fixed, as is done here. 

The SC results show a downturn for small to. This 
may be the beginning of an Anderson localization^* In 
the SC method, the phonons introduce a static (diag- 
onal) disorder which leads to an Anderson localization 
when the disorder becomes sufficiently strong. This is 
not expected to happen in the QMC calculation, since 
this calculation takes into account that the scattering is 
inelastic and phase information is lost J* We note, how- 
ever, that if such a downturn actually would occur in the 
QMC spectrum, it would be lost in the analytical contin- 
uation to real frequencies. This can be seen by comparing 
the SC and SCMEM curves, since the downturn in the 
SC results is lost in the SCMEM results, due to the MEM 
procedure. The reason is that the downturn happens on 
such a small energy scale that it is not detected by the 
MEM. 

Since there is a rather good agreement between the 
QMC and SC methods, we use the SC method in the 
rest of the paper. The reason is that the SC metod is 
simpler, and it is easier to interpret the results. 



RESULTS USING THE SEMICLASSICAL 
METHOD 



We have calculated the thermal conductivity of Nb us- 
ing the SC method. We used a cluster of 10 x 10 x 10 
atoms and averaged over 10 different distributions of 
thermally displaced atoms. The model described in Sec. 
ITT1 was used. The results are compared with experi- 
mental resultsi^iiLi 819 in Fig. Theory and experi- 
ment agree rather well. The figure also shows the curve 
k(T) = (47.5 + 0.013T) W/mK fitted to the experimen- 
tal results. This linear dependence on T is in qualitative 
agreement with the result expected from the discussion 
in Sec. I VII below. k(T) increases with T, which within 
the framework of the discussion in the introduction would 
correspond to the saturation of the apparent mean free 
path [see Eq. ©]. 

The thermal conductivity of A3C60 was calculated us- 
ing the SC method for a cluster of 8 x 8 x 8 molecules 
and averaging over 400 different distributions of ther- 
mally distorted molecules. The parameters were A = 0.5 
and u p h = 0.1 eV, where Wp^ = 0.1 eV is an approximate 
average renormalized phonon frequency for AsCeo- The 
results are shown in Fig. [21 The results are influenced by 
a downturn for small u> (see Fig. which depends on the 
amount of broadening used. For the temperatures shown 
in Fig.OK/csT 1 < 0.3 eV), however, this does not strongly 
influence the results. k(T) initially increases rapidly with 
T, but then reached a maximum followed by a rapid drop 
with T. This behavior differs qualitatively from the re- 
sults for Nb. This difference is discussed extensively in 
the next section. 
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FIG. 2: Thermal conductivity of Nb metal according to the 
SC method (full line) and from experiments ("Peletskiil"—-, 
"Moore"—, "Binkele"—, "Peletskii2"— ) together with the line 
k(T) = 47.5 + 0.013T adjusted to the experimental results. 
The calculated result was obtained for a 10 x 10 x 10 lattice 
using a Gaussian broadening with the FWHM of 0.06 eV. The 
results was averaged over 10 different random displacements 
of the atoms. The model of Nb described in Sec. ITIlwas used. 
The figure shows that theory and experiment agree rather 
well and that k(T) increases approximately linearly with T. 
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FIG. 3: Thermal conductivity of the Ceo model according to 
the SC method for A = 0.5 and ui p h =0.1 eV. 400 random 
configurations for a 8 x 8 x 8 cluster were generated, and a 
Gaussian broadening of 0.03 eV FWHM was used. 



Using the definition of }e in Eq. @ 

mE\l') = ~{e l -e l ,)Y,Ki{l\hi\l') 

i 

= +e l <)J2^(l\n<\l') (21) 

i 

= \fr + ei>m\l') 
It then follows that 

{i\h\i) = ^-Y L -mi\i')- (22) 

We now consider values of T which are sufficiently large 
to give a large thermal disorder. As a result there are ma- 
trix elements of the current operators between essentially 
all states. We therefore make the assumption that all the 
matrix elements of the particle current operator between 
different eigenstates of the Hamiltonian have the same 
value jav ■ This assumption has been checked extensively 
and found to lead to accurate results for a(ui) for large 
It is important to make this assumption for the par- 
ticle current operator, since it would obviously not be 
true for the heat current operator, as can be seen from 
Eq. due to the factor [(e + e )/2 - (j,)]. 

We first consider the case when T <C W . We can 
assume that the density of states per orbital and spin, 
N(e), is a constant N(n), since states close to e = fx 
mainly influence Im 7r nm (w) for small ui and T <C W. 
This gives 

ImTr^H = ^[NN d N{^) ]av ] 2 (23) 
x J dede'(^-ri 2 [f(e)-f(e')}5(hu; + e-e). 

We now consider the limit ui <C T. We also use the 
assumption T <§C W to extend the integrations in Eq. ll'MI) 
to infinity. To leading order in uj we then obtain 

Im^ 2) M = ~ 27r3( 3 fc * T)2 N[N d N{^)j av fn^. (24) 
Similar approximations lead to 

Im7r£ 2 >H = (25) 



VI. QUALITATIVE DISCUSSION IN 
SEMICLASSICAL FORMALISM 



In this section we use the semiclassical formalism to 
discuss qualitatively both T <C W and T ~ W, where 
W is the band width. For this purpose we first derive 
a relation between matrix elements of the electrical and 
heat current operators. In the spirit of the semiclassical 
approximation, we assume noninteracting electrons. 

Let \l) be a one-particle eigenstate with the energy 



and 



Im^M 



2tt 



NlNdN^jarftfcu. (26) 



From Eqs. (|12I13[) we then obtain 

«(0) = Tf_^B_Y 



a(Q)T 3 v e 
which is the Wiedemann-Franz law 



(27) 
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To obtain further understanding, we make some more 
explicit estimates. We first derive a sum rule for the 
particle current matrix elements. 



Hit' ,a\ 2 

ll'a 



NN d d 2 



(28) 



where a labels a coordinate, d is the separation of two 
atoms and (e 2 ) is the second moment of the density of 
states per site, orbital and spin. On the left hand side, 
the states \l) refer to extended eigenstates. In the middle 
expression, we have made a unitary transformation to lo- 
calized basis states and used Eq. Q, assuming nearest 
neighbor hopping only. Since the number of matrix ele- 
ments on the left hand side of Eq. (|2*5|l is 3(NN d ) 2 , we 
obtain 



ZNN d h z 



(29) 



Typically, (e 2 )N(n) 2 « 0.1 for a system close to half- 
filling. 7 Inserting this in Eqs. (|13l I24f> we obtain 



N d klT 

9 nn B 



k 2 T 



(30) 



where we have used that fi/d 2 = (4/3\/3)d, as appropri- 
ate for Nb metal with a bcc lattice. This gives the large 
T behavior k = 0.027T W/Km for Nb, where T is in K. 
This could be compared with the linear behavior fitted 
to experimental result 15 k = (47.5 + 0.013T) W/Km (see 
Sec. 0, where the linear term in T is about a factor of 
two smaller than our simple estimate. 

Comparing the result in Eq. H30|) with Eq. Q we find 
that we have to assume that 



I 



0.8N}, /3 d, 



(31) 



to obtain the estimated thermal conductivity. This is 
essentially the Ioffe-Regel condition. 

We now discuss the reasons for the saturation of the 
mean free path and the linear increase of n with T, using 
the semiclassical treatment of the phonons. At small T, 
there are important intraband transitions for small wave 
vectors q, which make contributions to Im n^ 22 ' 1 ^) for 
small u>. As T is increased, the increasing thermal dis- 
order leads to a strong violation of momentum conser- 
vation (in the electronic system) and many transitions 
that were strongly suppressed at small T become im- 
portant for large T. At the same time the transitions 
corresponding to small u are reduced. However, the sum 
rule in Eq. I|28(l shows that the sum over all particle cur- 
rent matrix elements squared is not reduced, as long as 
the hopping integrals are not reduced. Thus even if we 
assume that the particle current matrix elements corre- 
sponding to low energy transitions are not larger than 
the average of the matrix elements, there is still an ap- 
preciable (and saturating) weight of the transitions cor- 
responding to a small energy transfer in Eq. i|23|l . The 



prefactor [(e + e )/2 — /i] 2 in Eq. (|23|l furthermore gives 
increasingly large contributions as T is increased and the 
Fermi functions are broadened. This leads to the increase 
in k{T) with T and it corresponds to the increase in the 
specific heat per electron with T entering in the Sommer- 
feld theory leading to Eq. Q). 

We next consider the case when T ~ W . This is un- 
interesting for Nb, due to its large band width, but of 
relevance for A3C60 for which the band width is small. 
In A3C60 the phonons couple to the level energies. As a 
result the hopping integrals in Eq. (|28|l are not changed 
and j av is unchanged as T is increased. On the other 
hand, due to the fluctuations of the level energies, the 
band width grows with T as£ 



W(T) = W(T = 0)4/1 + c\ 



W{T = 0) : 



(32) 



where c is of the order 15. Thus there are substantial 
effects on the band width already for ksT ~ T^/10. To 
describe these effects we assume that the density of states 
can be written as 

N(e,T) = — 1 n( , £ ), (33) 

Vi + r/To Vi + r/To 

where n(e) is assumed to have no explicit T dependence 
and k B T rj 0.1 eV. We have 



. . 27re 2 N _ 
<j(uj) = — (N dJav ) 2 



de / de N(e,T)N(e ,T) 



x [/(e) -/(e )]8(hu + e-e). 



(34) 



In Rcf. I7J, a qualitative discussion of the electrical con- 
ductivity was given, neglecting finite T effects on the 
Fermi functions and emphasizing the difference between 
coupling to the level positions and hopping integrals. 
Thus the Fermi functions were replaced by ©-functions. 
If this is done in Eq. (|3*l|l and j av is assumed to be T 
independent, we obtain that cr(O) ~ 1/(1 + T/Tq) and 
the integral over cr{uS), entering the f-sum rule, goes as 
l/\/l + T/To, in agreement with the results in Ref. 0- 
For the treatment of the thermal conductivity, however, 
it is crucial to include the T dependence of the Fermi 
functions. The e integration can trivially be performed, 
leading to a factor /(e) — f(e + Tilo). Assuming that 
T ^S> W, we can approximate this factor as ?jw/(4fcsT). 
This gives 



cr(w) = 



--{N d j av f 



2iizy 1 + T/T Q 

x / dxn(x)n(x H =) 

; y/l + T/To' 

1 ha; 

g( n i=7=. )' 



(35) 



(T/ToVl + T/To y/l+T/To' 

i.e., for very large T, tr(0) decays as T~ 3 / 2 . Similar cal- 
culations for the thermal conductivity give 

k(u>) 



TThNy/T+T/% . 

-{N d j av ) 



2T 2 
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x / dxn(x)n(x + 



yfi + T / T o 

flLO 



2y/l + T /T 



(T/T ) 2 y/i + T/T c 



(36) 



As a result, the Wiedemann- Franz law is strongly vio- 
lated and k(0)/ct(0) ~ const instead of T for very large 
T. 

An important difference between T <C W and T ^> W 
is the behavior of the specific heat. For a fixed band 
width and a small T, the energy of noninteracting elec- 
trons increases as T 2 due to the smearing out of the Fermi 
function. This leads to an increase in the specific heat 
per electron, c v ~ T, as shown in Eq. and results in 
k/ct ~ T, For T > W, however, the electrons are already 
rather evenly distributed over the band width, and an 
additional increase of T does not increase the total en- 
ergy very much. The result is that c v then decreases with 
T. For such large T the semiclassical theory is invalid, 
and it is not possible to express k in terms of I and c v . 
Nevertheless, it is suggestive that k/ct ~ T is violated, 
although the proper power of T is not predicted by these 
arguments. 

To test the accuracy of these arguments, we define 



h(huj,T) 



(T/T 



y/l + T/T 



K(uy/T+T/lh), (37) 



and introduce an equivalent definition for g{hoj,T) in 
terms of a(u>). These functions are shown in Fig. 01 for 
different values of T. If the approximations introduced 
above were exact, the curves calculated for different val- 
ues of T would be identical. For the values of T shown in 
FigEl this requirement is quite well fulfilled for g(u,T) 
and less well fulfilled for h(ui,T). The main source of 
error is that the approximation for /(e) — f(e + hu>) in 
Eq. is only valid for fcgT ^> W, and Fig 0] shows re- 
sults for ksT > W/3. Nevertheless, varying T by a factor 
of almost ten only leads to a variation of the maximum 
values of h(u>,T) by less than a factor of two, providing 
support for the analysis above. 



VII. CONCLUSIONS 

We have studied the electronic part k of the ther- 
mal conductivity of metals, using two different method. 
Accurate results are obtained by using a determinantal 
Quantum Monte-Carlo (QMC) method together with a 
maximum entropy method (MEM). As a much simpler 
approach we use a method where the phonons but not 
electrons are treated semiclassically (SC). We applied 
these methods to a model of A3C60 (A= K, Rb) and 
showed that they give very similar results for n(u>) over 
most of the energy range if the renormalization of the 
phonon frequency is treated in the same way in both 
methods. For very small values of lo, however, the SC 
method gives a downturn. This is probably due to an 
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T=2.00 
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FIG. 4: Functions g(hw,T) and h(hco,T), defined in Eq. $7\ . 
as a function of w for different values of ksT — 0.25, 0.5, 1.0, 
1.5 and 2.0 eV calculated using the SC method. The parame- 
ters are A = 0.5 and cj p h = 0.1 eV. 400 random configurations 
for a 8 x 8 x 8 cluster were generated, and a Gaussian broad- 
ening of 0.03 eV FWHM was used. If the approximations 
behind Eqs. H35I36H were exact, the curves would fall on top 
of each other. 



incipient Anderson localization transition and a defect of 
the SC method neglecting that the scattering processes 
are inelastic. Applying the SC method to a model of Nb 
metal, we find that the results agree well with experi- 
ment. In particular, we find that k(T) increases with T, 
consistent with a saturation of the mean free path. In 
contrast, for A3C60 at very large T, we find a decrease of 
k(T) with T. The results are analyzed in the SC method. 
To discuss very large T we use an approximation where 
all matrix elements of the electric current operator are 
assumed to be the equal. We can then qualitatively re- 
produced the calculated results for Nb, and give a quan- 
tum mechanical understanding of why the apparent mean 
free path saturates. Within this framework it is also pos- 
sible to derive the Wiedemann-Franz law by assuming 
that T <§C W, where W is the band width. Due to the 
small band width of A3C60, we instead focus on the case 
T > W for this system. We then find that k(T) - T~ 3 / 2 
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indeed decreases with T and that the Wiedemann-Franz law is qualitatively violated. 
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